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Abstract.  Two  problems  of  Helmholtz's  equation  for  a  wave  incident  on  an  obstacle  are*  considered. 
For  the  first,  the  scattering  problem,  the  obstacle's  response  satisfies  Sommerfeld’s  outgoing  wave  radiation 
condition,  and  the  net  radiative  response  is  positive;  for  the  second,  the  response  satisfies  a  standing  wave 
condition  (an  appropriate  combination  of  outgoing  and  incoming  waves)  such  that  the  net  radiative  response 
is  zero.  The  essential  features  of  the  solutions  are  exhibited  in  terms  of  amplitude  functions  g  (the  usual 
scattering  amplitude)  and  g',  and  the  interrelation  of  the  functions  are  stressed  in  the  derivation  of  integral 
equations  g^g’^ (introduced  earlier  in  multiple  scattering  contexts).  The  scattering  amplitude  g  is  always 
complex,  but  the  simpler  function  g'  is  shown  to  be  imaginary  for  nonabsorbing  obstacles  having  inversion 
symmetry.  Long-wavelength  approximations  for  g'  may  be  obtained  from  potentiaUheory  and  perturbation 
procedures,  and  corresponding  approximations  for  g  then  follow  from  tW.-.  "S  ' 


!. 


Introduction.  We  consider  two  problems  for  Helmholtz’s  equations  correspond¬ 
ing  to  a  wave  incident  on  an  obstacle.  The  problems  differ  in  the  obstacle’s  response 
via  radiation;  for  the  first,  the  scattering  problem  (an  outward  radiating  obstacle),  the 
net  radiative  output  is  greater  than  zero;  for  the  second  the  net  radiative  output  is 
zero.  For  brevity,  the  obstacle  in  the  second  problem  is  labeled  radiationless  or 
nonscattering,  but  its  contribution  is  in  the  form  of  a  standing  wave  so  the  labels  refer 
to  net  energy  flow  via  radiation.  Section  1  defines  the  problems,  and  delineates 
analytical  and  physical  differences. 

The  solutions  of  the  problems  are  characterized  by  amplitude  functions  g  (the 
usual  scattering  amplitude)  and  g'.  Section  2  sketches  derivations  of  the  general 
theorems  for  g  with  emphasis  on  procedures  that  can  be  applied  for  g';  §  3  derives 
the  analogues  for  g',  and  §  4  applies  the  same  procedures  to  derive  integral  equations 
g[g']  and  their  inverses.  The  function  g'  is  simpler  than  g.  The  scattering  amplitude 
g  is  always  complex,  but  g'  is  shown  to  be  imaginary  for  lossless  (nonabsorbing) 
obstacles  having  inversion  symmetry.  Long-wavelength  approximations  for  g'  may  be 
obtained  from  potential  theory  solutions  and  perturbation  procedures,  and  correspond¬ 
ing  approximations  for  g  then  follow  from  g[g']. 

We  considered  g'[g]  =  g'[J;g]  initially  in  connection  with  multiple  scattering  of 
waves  by  periodic  arrays  of  obstacles  [1].  The  multiple  scattering  amplitude  (G)  of 
an  obstacle  in  the  array  was  related  to  its  response  (g)  in  isolation  by  a  functional 
equation  G  =  G[S;  g],  where  the  operator  5  equaled  a  discrete  infinite  sum  less  the 
analogous  integral.  It  proved  convenient  to  introduce  g'[J;g]  corresponding  to  g 
stripped  of  radiation  losses  and  to  express  G  functionally  in  terms  of  g'  and  a  modified 
operator  S'.  The  form  G[S';  g']  was  particularly  suitable  to  analyze  radiation  losses 
appropriate  for  the  array.  The  function  g  involves  radiative  losses  over  the  continuum 
of  real  directions,  but  G  involves  only  the  set  of  discrete  directions  corresponding  to 
the  propagating  modes  of  the  array  [1].  Working  with  g'  led  directly  to  simple  energy 
conserving  approximations  for  G,  and  to  rapidly  converging  representations.  The 
required  properties  of  g’  =  g'[f;  g]  were  deduced  from  known  properties  of  g. 
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Related  work  on  multiple  scattering  by  pair-correlated  distributions  of  obstacles 

[2]  led  to  functional  equations  G[g ]  involving  more  complicated  operators  than  S. 
To  reduce  these  to  tractable  forms  G[g']  for  nonsymmetrical  lossy  obstacles  required 
additional  properties  of  g'.  The  present  paper  derives  these  additional  properties  of 
g'  (as  well  as  those  obtained  earlier  from  g)  by  direct  considerations  of  the  radiationless 
ifbstacle  problem. 

The  following  includes  material  on  some  elementary  properties  of  the  scattering 
amplitude  g  and  the  associated  field  [3]— [8].  This  material  serves  to  distinguish  g  and 
g',  and  to  stress  that  the  nonscattering  problem  for  g'  does  not  correspond  to  a 
conventional  scattering  problem  subject  to  special  constraints  for  decreasing  the 
magnitude  of  g.  The  nonscattering  problem  we  consider  is  a  standing  wave  problem 
ancillary  to  the  scattering  problem;  similarly  for  the  associated  amplitude  g'  (labeled 
earlier  as  the  modified  or  transformed  amplitude)  and  g. 

1.  Statement  of  the  problems.  For  the  problems  at  hand  corresponding  to  a  plane 
wave  d»(r)  e~‘“"  of  angular  frequency  oj  incident  on  an  obstacle,  we  deal  with  solutions 
of  Helmholtz’s  equation  [3]-[8] 

(1)  (v2+fe2v  =  o,  *Ht| 

where  V2  =  V  •  V  is  Laplace’s  operator  and  k  is  the  propagation  parameter.  We  suppress 
the  time  factor  e~‘w\  but  all  waves  we  consider  have  implicit  period  2tt/id  in  t  and 
phase  velocity  w/fc;  the  plane  waves  have  explicit  period  27 r/k  (the  wavelength  A)  in 
the  space  coordinates. 

A  plane  wave  is  represented  by 

(2)  4>=e‘Vr,  k-kk,  r-ri 

with  directions  of  incidence  and  observation  given  by  unit  vectors  k  and  r.  In  three 
dimensions, 

(3)  r  =  r(0,  <p)  =  £cos  0  +  (xcos<p  +  y  sin  <p)sin  0,  k  =  r0  =  r(0o,  <Po), 

with  0  S  0  S  7r,  0  S  <p  <  27r ;  in  two  dimensions,  we  take  <p  =  0  and  0  §  0  <  2v ;  in  one, 
r  =  ±i,  and  0  =  0,  7r.  We  consider  three-,  two-,  and  one-dimensional  problems  in 
parallel  (corresponding  to  bounded  obstacles,  and  to  normal  incidence  on  cylinders 
and  slabs  respectively),  and  use  three-dimensional  terminology  for  all  cases.  If  three 
different  forms  of  a  function  or  parameter  arise,  they  are  sequenced  in  order  of 
decreasing  dimensionality.  The  center  of  the  obstacle’s  smallest  circumscribing  sphere 
(of  radius  a )  is  taken  as  the  phase  origin  (r  =  0). 

An  obstacle  (in  general,  a  discontinuity  in  the  parameters  associated  with  the 
medium  in  which  <t>  propagates)  is  specified  by  its  surface  9,  volume  V,  size  ka  relative 
to  wavelength,  and  by  conditions  on  in  5^  and  V.  An  obstacle  is  labeled  impenetrable 
if  ilr  =  0  in  Y  and  satisfies  on  9*  conditions  such  as:  either 

(4)  t£  =  0 

or 

(5)  A  •  =  dni//  =  0, 

with  n  as  the  outward  normal  into  the  external  volume  V,  or 

(6)  d„\}/  =  Z\l>, 

where  Z  is  constant  with  ImZSO.  (These  are  the  usual  Dirichlet,  Neumann,  and 
Robin  exterior  boundary  value  conditions  respectively.)  Physically,  e.g.,  for  small 
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amplitude  acoustics  with  ip  as  excess  pressure,  the  first  two  correspond  to  free  (pressure 
release)  and  rigid  surfaces  respectively,  and  the  third  is  a  generalization  in  terms  of 
the  impedance  Z. 

An  obstacle  is  labeled  penetrable  if  ip  satisfies  transition  conditions  on  y,  such  as 
(7  y)  <Ae  =  dn< h-  =Bdn4>i, 

where  E  and  I  stand  for  exterior  and  interior  (a  temporary  convenience),  and  i j/r  is 
a  nonsingular  solution  in  Y  of 

(7  V)  (V2  +  K2W=0,  K  ~  krjt  17 2  =  C/B. 

Here  the  obstacle’s  relative  parameters  C  and  B  are  constants  with  ImCgO  and 
ImfiSO,  and  17  is  the  corresponding  complex  relative  index  of  refraction.  In  acoustics, 
for  real  values,  C  is  the  obstacle’s  relative  compressibility  and  1  /B  its  relative  mass 
density;  complex  values  account  for  some  thermal  loss  effects,  compressive  viscosity, 
etc.  The  y  conditions  correspond  to  continuity  of  excess  pressure  (ij/)  and  normal 
excess  velocity  (proportional  to  n-  Vi£). 

Conditions  (4)— (7)  are  illustrative.  If  we  indicate  the  solutions  for  two  arbitrary 
directions  of  incidence  f  1,  r2  by  1 ,  (each  subject  to  the  same  conditions  on  y  and 
Y),  then  the  subsequent  development  applies  for  general  obstacles  for  which  Green’s 
surface  integral  for  <^1,  vanishes  on  y.  Introducing  a  symbolic  form  (and 
normalization  constants  for  subsequent  applications), 

(8)  {u,  u }  =  c  J  [u  (r')a„u  (t')-u  d„v  j  d^( r'),  c  = 

the  development  applies  for  all  cases  such  that  [7],  [8] 

(9)  «Az}  ==  0. 

By  inspection,  it  is  clear  that  (9)  covers  the  special  cases  (4)-(6).  To  show  that  (7)  is 
also  covered,  we  use  (7y)  and  then  apply  Gauss’  theorem  to  obtain  an  integral  over 
Y  which  vanishes  by  (7  Y). 

The  conditions  on  the  nonvanishing  imaginary  parts  of  Z,  C,  B  correspond  to 
energy  absorption  by  the  obstacles  (lossy  obstacles).  Changing  the  signs  of  the 
imaginary  parts  corresponds  to  energy  production  (gainy  obstacles,  or  negative-lossy 
obstacles).  If  satisfies  (4)  or  (5),  or  (6)  for  real  Z,  or  (7)  for  real  B  and  C,  we  label 
it  as  lossless.  More  generally,  a  lossless  obstacle  subject  to  (9)  also  satisfies  [7],  [8] 

(10)  c^z}  =  0, 

where  t/r*  is  the  complex  conjugate  of  t p.  The  statement  after  (9)  applies  equally  for  (10). 
In  general,  (10)  is  not  satisfied.  If  we  normalize  the  incident  energy  flux  vector  as 

(11)  Re  (4>  *V(j>/ik )  =  k, 

with  unit  energy  density  in  unit  area  normal  to  k,  then  for  lossy  obstacles  the  net 
energy  flow  into  y  (the  absorbed  energy)  is  given  by 

(12)  -J  Re ((^*V^/i^)-fi dy  = -2<ro{^*, cr0=l/icfc, 

where  crA  >  0  is  the  absorption  cross  section.  Were  <rA  <  0,  we  would  relabel  it  -<tp 
with  <7P>0,  and  cancel  the  minus  signs  to  discuss  energy  production  (negative- 
absorption). 
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The  two  classes  of  problems  we  consider  are  determined  by  different  conditions 
for  i/»  -  <(>  on  the  surface 

(13)  5  =  lim  S(r). 

r— co 

We  label  S  as  the  spherical  surface  at  infinity,  the  volume  V  external  to  as  the 
shell  volume,  and  suppress  limit  operations  in  general. 

One  problem,  the  scattering  or  radiating  obstacle,  is  well-posed  if  we  supplement 
(1)  and  (9)  by  the  Magnus-Wilcox  form  of  Sommerfeld’s  radiation  condition  [3],  [4] 
fori/f-<£  =  m(  r), 

(14)  /  =  J  \drU-iku\2dS  =  0. 

This  condition  leads  to  [4] 

(15)  \\u\2dS  =  b,  Jm  2dS**b' 

with  b  and  b'  finite,  and  to  [4] 

(16)  |  Re  (u*dru/ik)  dS  =  J  Re  (u*Vu/ik)‘hdy  =  2(ro{u*,  u}^<rs>0, 

where  as  is  the  scattering  cross  section.  Equation  (16)  states  there  is  net  energy  outflow 
from  the  obstacle  [3],  [4].  The  associated  wave  ue~,U1'  will  be  shown  to  have  phase 
kr  —  tot  for  r~ 00  corresponding  to  a  wave  u*  e~,a"  outgoing  from  r  =  0.  The  sum 
( r  =  (Ta+<ts,  the  total  energy  cross  section,  will  be  shown  proportional  to  the  interfer¬ 
ence  terms  of  </>  and  u  in  the  forward  direction  r  =  k.  Thus  cr  is  a  measure  of  the 
energy  the  obstacle  derives  by  interference  of  <f>  and  u  and  dissipates  via  absorption 
(conversion  to  heat,  to  radiation  with  different  w,  etc.)  or  diverts  as  scattered  radiation 
over  the  continuum  of  directions  f. 

Because  crs>0  in  (16),  the  energy  is  outflowing  from  the  obstacle  which  thus 
represents  a  secondary  source  of  radiation.  Were  we  to  change  the  sign  in  (14),  we 
would  obtain  as  <  0.  and  energy  would  flow  into  the  obstacle  corresponding  to  a  sink 
of  radiation.  The  associated  wave  u~  e~“"  would  have  phase  -kr-ut  for  r~ 00 
corresponding  to  a  wave  incoming  to  r  =  0. 

The  second  problem  we  consider,  for  d/  =<//'  =  d>  +  u'  subject  to  (1)  and  (9),  is 
that  of  the  radiationless  obstacle 

(17)  J  Re  (u'*dM'/ikiJ‘*0, 

(18)  u'  =  i(u  +  -u“),  I^llB^iku^dS-O, 

where  u'  represents  a  standing  wave,  an  appropriate  combination  of  an  outgoing  and 
incoming  wave  corresponding  to  zero  net  radiation.  The  S  integral  in  (18)  gives  (15) 
for  n*,  and  erf  2  0  for  (16).  In  order  to  facilitate  discussion  of  u\  we  first  consider 
the  radiating  obstacle  specified  by  u  subject  to  (14).  The  development  for  u  covers 
aspects  common  to  both  u  and  u'  with  emphasis  on  procedures  that  apply  for  both, 
but  some  elementary  properties  of  u  and  its  associated  scattering  amplitude  g  are 
included  to  delineate  differences  between  u  and  u'. 
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2.  The  scattering  obstacle;  the  amplitude  g.  We  apply  Green’s  theorem  to  u  and 

(19)  afo(fc|r-r'|)  =  af0(p)  =  /,<%),  Z/'11  (p),  *“>, 

where  the  appropriate  zeroth  order  Hankel  function  of  the  first  kind,  is  the 
normalized  free-space  Green’s  function.  We  consider  a  fixed  observation  point  r  in 
the  shell  V (r'),  exclude  the  singularity  of  X0  by  a  limiting  sphere  with  radius  |r  -  r'|  -» 0, 
and  obtain  u  as  the  difference  of  Green’s  surface  integrals  over  if  and  S 

(20)  u  (r)  =  {%o (k  |r-r'|),  u(r')}-c  j  (X0dnu  -  u  d„XQ)  dS. 

We  proceed  essentially  as  in  [4]. 

The  integral  over  S  can  be  written  in  terms  of 

/i  =  J  X0 (d„u  - iku )  dS,  I2  =  |  Xdu/r')  dS,  %0 - X(kr') 

where  |$f|<x(r')  <m-l)/2  and  dSoc(r')m_1  for  m=3,2,  1;  /2  does  not  arise  for  the 
one-dimensional  case,  and  the  limit  for  r'-»oo  is  implicit.  Using  Schwarz’s  inequality, 
and  comparing  |/i|2  with  /  of  (14)  and  |/2|2  with  b  of  (15),  we  obtain  f(-)rfS  =  0  in 

(20) .  Thus 

(21 )  u  (r)  =  {%0{k  |r  -  r'|),  u  (r')J  =  {%Q,  u }, 

with  the  brace  operation  over  if{ r')  as  in  (8)  or  over  any  surface  that  separates  if 
from  r,  is  a  radiative  solution  of  Helmholtz’s  equation  [3],  [4],  For  r  —  oo, 

(22)  u(r)-~9f(kr)g(i,  k), 

(23)  g(f,  k) ={e u (r')}  =  u)  =  {^f,  u), 

with  Sf  as  the  asymptotic  form  of  #?0,  and  g(r,  k)  as  the  normalized  scattering  amplitude 
[7],  [8].  The  function  differs  from  2C0  only  in  two  dimensions,  i.e., 

e'°  /  7  \  U2 

(24)  X{p)  =  — ,  (  — )  eip~iv,\  eip. 

ip  \irpJ 

In  one-dimensional  problems  we  use  p/k  =  |z|  =  ±z  and  g(±*,z).  For  all  cases  we 
write  the  phase  of  the  associated  wave  as  kr  —  a>t  corresponding  to  a  wave  outgoing 
from  r  =  0. 

The  scattering  amplitude  g  not  only  determines  u  for  r  ~  oo  (i.e.,  in  the  far-field 
kr  » 1,  r  »a)  but  for  all  r>a,  and  at  least  for  some  :ia,  Substituting  Noether’s  [6] 
and  Sommerfeld’s  [5]  complex  spectral  representation  for  X0  into  (21)  and  using  the 
definition  (23)  for  g,  we  obtain  [8],  [7]  a  corresponding  spectral  form  for  u.  Thus,  at 
least  for  r>a  for  all  ?,  and  for  r  greater  than  the  obstacle’s  projection  on  r, 

(25)  M(r)=fe'^-,g(f„t). 

Jc 

In  three  dimensions  [8],  (c  =  t(6c,  <pc )  and  jc  =  (l/2rr)  J  dCL(Bc,  <pc)  with  contours  [6]  as 
for  ho  * ;  in  two  [7],  fc  =  ?(0C)  and  jc  =  (l/w)  f  ddc  with  contour  [5]  as  for  Ho' ;  in  one, 
fc  =  ±i  and  |c  selects  the  sign  corresponding  to  z  =  ±|z|.  (The  associated  integral 
operator  for  /0  =  Re  X0  is  the  mean  (M)  over  real  directions  of  observation.)  The 
brace  operator  in  (23)  specifies  g  in  terms  of  u,  and  the  fc  operator  in  (25)  is  an  inverse 
for  u  in  terms  of  g.  Other  inverses,  and  complete  convergent  and  asymptotic  expansions 
of  u  for  large  kr  (in  terms  of  g  and  its  derivatives  with  respect  to  angles)  are  discussed 
in  the  literature  [7}— £10]. 
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From  the  general  obstacle  condition  (9),  in  terms  of  ip,  =  <pi  +  uh  we  have 

(26)  {ip  i ,  tp2)  =  {<P  l  +  u  i ,  <p2  +  m2}  =  {<£  i .  m2}  +  {« 1 .  <£2}  =  0, 

where  {<^1,  <£2}  vanished  by  Green’s  theorem  and  (V2  +  k2)<pi  -  Oin  V.  The  term  {u  1 ,  u2} 
vanished  by  using  Green’s  theorem  to  convert  the  PP  integral  to  V  plus  S  integrals; 
the  V  integral  is  zero  from  (1)  for  u,  and  the  S  integral  is  zero  because  the  asymptotic 
form  u  of  (22)  is  appropriate.  From  (26),  we  have  {<p u  u2)  =  {<f>2,  Mi};  by  the  definition 
of  g  as  in  (23)  this  corresponds  to 

(27)  g(-fi,f2)  =  g(-f2,?i),  g(?,£)  =  g(-£,-?), 

i.e.,  to  the  usual  reciprocity  theorem. 

Proceeding  similarly  for  {ipf,  ip2},  we  have 

{l Pi,  lpi)-{<P\,  «2}  +  {m*,  <£2}  +  {m*,  «2} 

(28) 

*  g(?  1.  ?2)  +  g*(?2,  ?  1)  +  C  j  (u*dnu2-  U2dnu  1 )  dS, 

where  we  used  {«*,  <p2)  =  ~{<p2,  uf}  =  {<p*,  Ui}*.  In  the  integral  over  S,  we  substitute 
the  asymptotic  u  of  (22)  to  obtain 

(29)  {ipf,  ip2}  =  g(ru  r2)  +  g*(r2,  ii)  +  2J(g*(r,  fi)g(r,  f2), 
where  M  is  the  mean  over  all  values  of  r.  Explicitly 

(30)  =  f  da(6,<p)  =  ~  f  d<p  f  ddsmS,  M2  =  ^~  [  de, 

4  tt  J  4  7T  Jo  Jo  27t  Jo 

and  M 1  is  one-half  the  forward  and  back  values. 

For  lossless  obstacles,  (10)  requires  {<p*,  tp2}  =  0  and  (29)  reduces  to 

(31)  -g(r„  ?2)-g*(r2,  fi)  =  2Mg*(r,  r,)g(f,  r2). 

In  particular,  if  the  obstacle  has  inversion  symmetry  (so  that  ron^  implies  -r  on 
y),  then  from  (27)  we  have  g(iu  r2)  =  g(?2,  fi),  and  from  (31), 

(32)  -Re  g(f„  f2)  =  Atg*Q,  f  ,)g(f,  f2); 

for  such  symmetry,  the  right  side  is  real.  If  |Im  g) » jRe  g|,  then  Re  g  is  of  order  (Im  gf 
and  (32)  may  be  used1112  to  construct  Re  g  from  approximations  of  Im  g. 

From  (31)  for  ?i  =  ?2  =  k,  we  obtain  the  forward  scattering  theorem  with  no 
requirement  of  symmetry, 

(33)  -Re  g(fi,  ft)  =^|g(f,  fc)|2  =  Mu*,  m}  =  crs/<r0 

in  terms  of  crs  of  (16).  This  relation  corresponds  to  the  energy  theorem  for  lossless 
scatterers:  the  energy  derived  by  the  obstacle  via  interference  of  the  incident  wave  <p 
with  the  scattered  wave  u  in  the  forward  direction  fi  is  radiated  over  all  directions  i. 
The  theorem  also  states  that  if  Reg(fi,  £)  =  0  (if  there  is  no  interference  of  <p  and 
n),  then  g(f,  fi)  =  0  for  all  ?  (then  there  is  no  scattering):  Re  g(fi,  fi)  vanishes  only  if 
^f|g(f,  fy)2  vanishes,  and  this  requires  Re2g(f,  £)  +  Im2g(f,  £)  =  0,  from  which 
Reg(f,  fi)  =  Img(f,  fi)  =  0.  Since  fi  is  arbitrary,  Reg(fi,  fc)  =  0  implies  g(f,  £)  =  0  for 
all  ?. 
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For  lossy  obstacles,  from  (29)  in  terms  of  (12)  and  (16), 


(34) 


-Re  g(k,  ft)  =  -2{<A*,  0}+j{w*>  «}  =  (o-A  +  <7s)/o- o, 


4rr  4 

0-0  =  F’?2’ 


which  states  that  the  energy  derived  by  the  obstacle  via  interference  of  <f>  with  u  is 
balanced  by  the  energy  absorbed  and  reradiated.  Because  aA  and  <rs  are  positive, 
both  must  vanish  in  order  for  Re  g(k,  k)  to  vanish;  if  there  is  no  interference,  there 
is  neither  absorption  nor  scattering  (no  obstacle).  On  the  other  hand  for  a  gainy 
obstacle,  such  that  excitation  by  4>  causes  it  to  produce  energy  -crA  =crP>  0,  then 

(35)  — <r0  Re  g(k,  k)  =  -aP  +  as 

may  vanish:  if  crP  =  trs,  no  energy  is  contributed  via  interference  of  <t>  and  u  ;  the 
energy  trP  produced  by  the  obstacle  balances  the  energy  it  radiates.  Although 
Reg(k,  k)  may  vanish  for  o-s  #  0  in  (35),  the  converse  is  ruled  out  (except  for  the 
trivial  case);  the  vanishing  of  as  =o-<y#|g(r,  k)j2  implies  Reg(k,  k)  =  0,  so  that  o> 
would  also  vanish  (no  obstacle).  If  aP>as,  then  Reg(k,  k)>0  and  the  constructive 
interference  between  <f>  and  u  is  a  result  of  the  energy  produced  by  the  obstacle. 

For  the  general  case  specified  by  (28),  the  analogue  of  (34)  is 

-g(f i»  h)  ~ g*(r2,  f i)  =  -{i/'f,  ifri}  +  {«  f,  u2}  =  2(o- a1  +  a\ 1  )/a0, 

^  as 1  =  a&Mg*(x,  rOg(r,  r2). 


If  the  obstacle  has  inversion  symmetry,  the  left  side  reduces  to  -2  Re  g(f1(  r2)  as  in 
(32),  and  a2A  and  al 1  are  real. 

3.  The  nonscattering  obstacle;  the  amplitude  g'.  We  consider  the  problem  for 
<A’  =  d>  +  u'  defined  by  (1),  (2),  (9),  (17)  and  (18).  The  problem  for  i //'  differs  from  the 
scattering  problem  for  i/»  in  that  the  radiation  condition  (14)  designed  to  correspond 
to  net  outgoing  radiation  (16)  is  replaced  by  the  nonradiating  (17)  standing  wave 
(18)  conditions.  Relabeling  the  outgoing  Green’s  function  (19)  and  introducing  its 
complex  conjugates  as  the  incoming  one, 

QtS>  _  at/?*  _  qu>~  _  q^2) 

/ )  <ftO  —  erto  —  £7t  0  t  (Tt  0  —  at  o  — 

we  see  that  the  standing  wave  Green’s  function 

(38)  te -^o)  =  i'Im^0  =  Uf. o 
has  the  same  singularity  as  #f0. 

Applying  Green’s  theorem  to  jjV0  and  and  proceeding  as  for  (20),  we  obtain 
u '(r)  =  {uVoik |r - r'(),  u V)} -(/,,+  hr ~ hr -In ), 

(39)  u'  =  k(u  +  -u~)  =  k(un)-ua\ 

4 /„  =  c  J  (*U ?dnu  u)  -  u  U)dn&o)  dS. 

The  integrals  I\  i  and  /22  can  be  expressed  in  terms  of 

IX  =  |  xs  (dnu *  T iku ±)dS,  It  =  |  #0  (u ±/r')  dS ; 
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comparing  |/f  |2  with  /*  of  (18),  and  |/|  \2  and  b  of  (15),  shows  that  In  and  I22  vanish 
essentially  as  discussed  for  (20).  The  integrals  / 12  and  I2\  involve  analogues  of  I2  in 
terms  of  u  *  which  also  vanish  by  comparison  with  b,  so  that  the  key  integrals  are 

|  %o  (dnu *  ±iku*)dS  =  J±  ± i2*  |  SHo u *  dS, 

where  J±,  differing  from  It  only  in  the  interchange  of  %+  and  #f~,  vanishes  by  (18). 
Consequently,  to  eliminate  f(  • )  dS  of  (39)  we  require 

(40)  R  =  j{3tfiu*-XZu-)dS  =  0 

which  constrains  the  amplitudes  of  u  +  and  u  Thus  subject  to  (40),  we  have 

(41)  u '(r)  =  {LN"o(k |r - r'|),  n'(r')}  =  l{^,  u'}-\{% 0,  «'} 

For  r  ~  00, 

(42)  n'(r)~^+(fcr)g'(f,k)-^-(ifcr)g'(-r,k), 

(43)  g'(f,k)  =  {e-ifc,v,H'} 

where  $?+  =  VC  of  (24),  =  9f*,  and  g'  is  the  amplitude  for  the  radiationless  problem. 

For  all  cases,  u' e~"“'  consists  of  outgoing  plus  incoming  waves  with  corresponding 
phases  kr  -  oit  and  -kr  -  tot. 

To  ascertain  the  import  of  (40)  we  consider  u*  further.  From  In  =  /2 2  =  0,  and 
from  the  signs  associated  with  $fo  (p)/i  for  p  -»0,  we  have 

(44)  «lr)=±{^o,K*},  u*~±X±{e*ik,t\u±}  =  ±X*g*(±i), 

and  consequently 

(45)  «'(r)  =  i(„  +  -  « 1  ~  g+(r)  +  Wg  1-f), 

where  k  was  suppressed  as  unessential.  Alternatively,  from  (41)  with  u  in  the  braces 
replaced  by  |(u+-u“), 

(46) 

«'  ~  i*Tglf)  -  gif)]  -  [gl-f)  ~  gl-f)]. 

From  (45)  and  (46), 

(47)  glf)  =  5[glf)-glf)],  -gl~f)  =  2[gl-f)- gl-f)], 
so  that  the  amplitudes  are  related  by 

(48)  g+(f)  =  -gif) 
for  all  f.  From  (42),  (45)  and  (48),  we  also  have 

(49)  g'(f,  k)  =  g+(f),  g'(-f,fe)  = -gl-f)  =  gl-f) 

Substituting  (44)  into  the  constraint  R  =  0  of  (40),  we  see  that  it  corresponds  to 

(50)  R  =  j  ^T^r[^glf)+4,_(g-(-f)]r/5  =  0, 
or  equivalently,  to  the  requirement  that 

(51)  fc)  -  <*->g'(-f,  fc)]  =  0. 
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That  the  mean  value  in  (51)  vanishes  follows  from  the  general  form 

(52)  -/(-?)]  =  o, 

i.e.,  r  is  a  dummy  and  the  operator  M  means  integrate  over  all  r.  Thus  the  constraint 
R  =  0  is  satisfied  and  the  procedure  is  fully  consistent. 

We  now  apply  (41)— (43)  to  parallel  the  development  (26)-(36).  Because  <A'  = 
<£  +  u'  satisfies  the  general  obstacle  condition  (9),  we  have 

(53)  {*li\,<l>2)  —  {4>i>U2}  +  {u'u  <t>-iS  +  {u'u  u'2)s  =  0, 

where  {4>i,  <fr2)  vanished,  and  {u\,  u2)  was  converted  to  an  S  integral  via  the  discussion 
for  (26).  On  S,  we  use  (42)  to  reduce  the  integral  to 

(54)  ^f[g'(f,  ?i)g'(-r,  f2)  -  g'(-r,  ?i)g'(r,  r2)]  =  0, 

which  vanishes  by  (52).  Thus  (53)  reduces  to  {<^1,  m2}  +  {m'i,<£2}  =  0,  and  proceeding 
as  for  (27)  in  terms  of  (43),  we  obtain 

(55)  g'(-?i,f2)  =  g'(-*2,fi),  g'(r,  k)  =  g'(-k,  -r). 


i.e.,  the  same  reciprocity  theorem  as  for  g. 

Proceeding  similarly  for  {<l/\  *,  02|,  we  construct  the  analogue  of  (28), 

(56)  {iA'*,^2}  =  g'(fi,  r2)  +  g'*(r2,  r, )  +  {«'*,  Mis¬ 


using  the  asymptotic  u'  of  (42)  on  S, 


(57) 


{ui*,  u’2}s  =  f,  r,)g'(r,  f2)  -g'*(-f,  f,)g'(-f,  h))  =  0, 


which  vanishes  by  (52).  The  special  case  of  (57)  for  f,  =  r2  =  k  corresponds  to  the 
analogue  of  the  scattering  cross  section: 


(58) 


O-i  =  CToU  |g'(f,  k)|2  -  |g'(-f,  k)|2]  =  0. 


Thus,  the  obstacle  is  radiationless  in  the  sense  required:  there  is  no  net  outflow  or 
inflow  of  energy  in  the  form  of  radiation  over  the  continuum  of  directions  r.  Con¬ 
sequently  (56)  reduces  to 

(59)  -g'(tu  r2)  -  g'*(t2 ,  f ,)  =  -{<*;*,  =  2<ta(?2,  f l)/cr0, 

i.e.,  to  the  analogue  of  (36)  for  no  scattering  losses.  In  particular,  for  r  1  =  r2  =  k, 

(60)  -Re  g'(k,  k)  =  -\W*,  <£'}  =  o-^/o-o, 

where  a\  is  the  absorption  cross  section  of  the  radiationless  obstacle.  The  energy 
theorem  (60)  states  that  the  energy  the  obstacle  derives  via  interference  of  <t>  and  u' 
in  the  forward  direction  is  dissipated  by  absorption.  The  presence  of  the  interference 
term  qualifies  the  label  radiationless  introduced  for  brevity:  the  label  radiationless  is 
defined  by  (17)  and  (58). 

If  the  obstacle  is  lossless,  then  (59)  reduces  to 


(61) 


g'(fl,*2)  +  g'*tf2,?l)  =  0. 


If  the  obstacle  has  inversion  symmetry, 

(62)  Re  g'(?i,  f2)  =  0 


and  even  if  it  doesn’t,  the  forward  value  of  (61)  satisfies 


(63) 


Re  g'(£,  fc)  =  0. 
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There  is  no  interference  of  <f>  and  u'  for  this  case;  the  obstacle  neither  scatters  nor 
absorbs,  so  no  energy  is  called  for. 

4.  Interrelations  of  g  and  g\  The  general  obstacle  condition  defined  in  (9)  was 
applied  as  an  operator  to  two  solutions  tl/u  >p2  for  the  scattering  obstacle  in  (26)  and 
to  two  solutions  i/ri,</r2  for  the  nonscattering  obstacle  in  (53).  This  led  to  the  same 
reciprocity  relation,  (27)  for  g  and  (55)  for  g',  for  the  values  of  the  amplitude  with 
arguments  involving  two  arbitrary  directions  of  incidence  ?i,  r2.  We  now  apply  the 
operator  to  ip\  and  <^2>  where  </»  and  <A'  satisfy  the  same  conditions  on  &  and  V.  This 
leads  to  an  integral  equation  for  g  in  terms  of  g’,  as  well  as  to  the  inverse. 

As  discussed  for  (26)  and  (53), 

(64)  {l/'l,lA2}={<f>l+Kl,^2  +  «2}  =  {<^l.  H  2  }  {<$ 2 >  Kl}  +  0<1,  «2}s  =0. 

Using  the  definitions  of  g  and  g'  as  in  (23)  and  (43),  and  the  asymptotic  u  and  u'  as 
in  (22)  and  (42),  we  obtain 

g(-?2,  fi)  =  g'(-?i,  r2)  +Mg(  r,  fi)g'(-f,  r2) 

=  g'(-?2,  f ,)  +  Jtg'(-i2,  r)g(f,  f ,), 

where  the  second  version  follows  from  the  reciprocity  theorem  for  g'.  Relabeling  -r2 
as  r2,  we  have 

(66)  g(f2,  f ,)  =  g'(r2,  f,)  +Mg'(r2,  r)g(r,  r,) 

which  applies  as  an  integral  equation  for  g  in  terms  of  g',  or  for  the  converse.  This 
form  and  a  related  one  arose  earlier  [1]  in  reducing  functional  equations  for  multiple¬ 
scattering  problems  which  involved  an  operator  consisting  of  M  plus  additional 
operations.  The  related  form  follows  from  (66)  by  using  (52)  as  Mf{r)  -  Mf(-i)  and 
then  applying  the  reciprocity  theorems.  We  construct  ^fg'(r2, -f)g(-f,  fi)  = 
Mg'(i,  -r2)g(-fi,  f)  =  g(-?i,  -r2)-g'(-fi,  -r2),  and  then  relabel  -?,  as  f,  to  obtain 

(67)  g (f i ,  r2)  =  g'(fi,  f2)  +J(g (iu  f)g'(r,  f2). 

Since  fi  and  r2  are  arbitrary,  the  only  substantive  difference  in  (66)  and  (67)  is  in  the 
location  of  f.  Form  (67)  follows  trivially  from  (66)  if  the  obstacle  has  inversion 
symmetry,  but  the  present  procedure  shows  it  holds  for  all  shapes. 

Similarly  we  apply  the  generalized  energy-type  operation  considered  in  (28)  for 
■A  and  (56)  for  </r'  to  the  mixed  set  <A'*,  < fa  to  obtain 

{<a  ^ ? ,  u2} + {u  <f>2) + {u  u2}s 

(68) 

= g(f  1,  *2) + g'*(?2,  fo +Mg'*d,  fog(f,  f2). 

For  lossless  obstacles, 

(69)  g(f  1,  r2)  =  -g'*(f2,  r  1)  -^fg'*(r,  ri)g(f,  f2), 

as  follows  directly  from  (66)  by  using  (61)  for  lossless  g'.  From  (69)  and  (61)  (or  from 
(66),  (67),  and  (61)),  we  have 

-*(».,  *2>-g*(fa,  ft)  =  M[g'*d,  f,)g(f,  f2)  +  g'(f,  f2)g*(f,  f,)] 

(70) 

=  ^[g'*(f,  f1)g(f,  f2)  +  g'(ft,  f)g*(f2,  f)]. 

If  the  obstacle  has  inversion  symmetry,  then 

(71)  -Re  g(f,,  f2)  =  Re  Mg'*(t,  f,)g(f,  f2)  Im  g'(f,  f.)  Im  g(f,  f2), 
where  g'  =  i  Img'. 
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For  nonsyumetrical  obstacles,  from  (66)  and  (67)  in  briefer  notation 

(72)  gl2  =  gl2  +  ^glrgr2.  gfl  =  g  21  (g2rgr  1  )*■ 

Using  (66)  for  lossless  obstacles,  i.e.,  g2*  =  ~gn,  we  convert  g*i  to 

(73)  g?i  =-g[2~J(gugZ, 
which  together  with  gi2  of  (72)  leads  to 

(74)  gl2  +  g*l  ="^glr(gr2~  g*r),  gl2  ~  g  21  =  2g  1.2  +Mg  [,  (gr2  +  g  *>■). 

Writing 

_  gl2  =  Si2  +  A12,  g21  —  Si2~Ai2,  5i2  =  5(gl2  +  g2l).  A  12  =  |(gl2  ~  g2l). 

(75) 

P12  =  gi2  +  g*i  =  2  Re  Si2+2i  Im  An.  W12  =  gi2  ~  g*  1  =  2i  Im  S12  +  2  Re  An, 
We  recast  (74)  as 


(76) 


W12  —  2gn  +Mrg'lrpr2>  P 12  —  M,g\rmr2- 


Thus,  introducing  M.p,  we  uncouple  m  and  p  by 

fti  12  2gn  JrMrguMpgrpmpl, 

(77)  (  (  t  t 

Pi 2  2-7/ |  rgr  2  "^■/^i-glr4^pgppPp2, 

which  supplement  the  original  integral  equation  (66)  when  the  obstacle  is  lossless. 

5.  Supplementary  considerations  and  applications.  The  form  (67)  arose  earlier 
for  the  two-dimensional  problem  of  scattering  by  a  grating  of  equal  spaced  parallel 
cylinders  [1].  Using  a  symbolic  procedure  and  the  energy  theorem  (32)  for  lossless  g 
having  inversion  symmetry,  it  was  shown  that  Re  g \2  —  0  of  (62)  was  satisfied.  Writing 
(67)  for  the  two-dimensional  case  (with  ?>,  ?2  replaced  by  6 1,  02  with  no  function¬ 
relevant  connotation)  as 


(78) 


g(#l,  02)  —  g'(01,  #2)  +« —  [ 

2rr  Jo 


2n 


g(0u0)g'(0,02)d0 


2ir 


|  d0[2n8(0 


~0i)  +  g(0i  ,0)]g'(0,02) 


provides  the  starting  point  [1,  (55)]  of  the  earlier  development.  The  corresponding 
symbolic  version  in  [1,  (55a)], 


(79)  Pg'  =  g,  g'  =  P“’  g, 
was  used  to  construct  [1,  (57)] 

(80)  g'  =  p-'g  =  |p-fp*g  =  |  P'T  [g(0,  02)  +  ~  j  g*(0,  0')g(0\  02)  dtf] 
Applying  the  two-dimensional  case  case  of  (32)  led  to 

(81)  Reg'(0i,  02)  =  0,  g'(eu02)  =  \P~1\2ilmg(6,02) 


where  the  first  equality  corresponds  to  (62);  for  the  second,  we  invert  via  P*(Pg')  =  P*g 
to  reconstruct  (32).  See  [1]  for  applications  of  (78)  and  (81)  in  reducing  the  multiple 
scattering  amplitude. 
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The  integral  equation  (78)  was  also  applied  in  [1,  (104)]  to  circular  cylinders  with 
amplitudes 

(82)  g(0,a)  =  l  g'(6,  a)  =  £  a'ne""e-a) 

rt  —  —  ac 

to  obtain 

(83)  a'„=an-a'nan,  a'n  =  a„/(\ +a„) 


in  terms  of  the  known  scattering  coefficients  a„.  Using  g  in  the  energy  theorem  (32) 
for  lossless  scatterers  led  to  -Re  an  =  |«„|2,  and  consequently  to 

(84)  a'n  =  i  Im  <z„/|l  +  an\2 

as  discussed  for  [1,  (104)]  and  exhibited  for  the  boundary  conditions  (4)  and  (5). 

The  coefficients  an  and  a'n  corresponding  essentially  to  all  conditions  (4)-(7)  were 
considered  before  [11]  in  the  context  of  scattering  by  random  distributions  of  radially 
symmetric  obstacles.  Thus,  all  cases  are  covered  by 


(85) 


an 


Rnjn  ~  fn 

n  q//>t  ’ 

^  n 


=J’n(ka ),  +  X" 


where  $  =  /,  /  and  31?  =  A,  H,  and  the  prime  indicates  differentiation  with  respect  to 
argument.  (See  [11]  for  one-dimensional  analogues  for  31?  for  n  =0, 1.)  for  the 
boundary  conditions  (4),  (5),  (6)  we  let  R  -  oo,  0,  Z/k  respectively;  for  the  transition 
conditions  (7), 

(86)  Rn~zrAKa)l$n(Ka),  (~Br,=(CB)u2 

where  £  represents  an  impedance. 

From  (83),  we  have 


(87) 


a'n=i(RJn-'?'n)/(Rn^n-Srn)^i(Rn?n-fn)/Dn 


where  and  Jfn  are  real  functions  of  ka,  and  Rn  is  real  if  the  parameters  Z  or  B 
and  K  are  real  (lossless  obstacles).  In  general 

(88)  Re  ai  -  W  Im  RJ\Dn\z  =  -H^|Im  Rn|/|£>n|2,  W  =  l/(kaf,2/7rka, 
where  W  =  $JC'n  is  the  Wronskian.  Inverting  the  development  of  (84), 


(89) 


an  =; 


^  rt 


a'n~\a'n\2 

r. 


1  -a'n  |l-fl'| 
which  for  lossless  obstacles  yields 

(90)  -Re  a„  = 

More  generally 


«n+|«n|2  =  pj - VjT. 

1 1  —  a  n  I 


la;12 


i  +  K 


T7i~\an\  ■ 


(91) 


r>-  .  i  _  |2  a'n 

Re  an  )an  |  —  r:  rfj, 

|1  -a„\ 


,  ,2_  |a'|2 

\an\  = 


with  -Re  a'n  as  in  (88)  isolates  the  effects  of  absorption. 
For  all  such  problems  we  have  the  forms 


=  —  jg  a'n 

n  Hn  Jn+fMn  =  l-a'’ 


(92) 


a, 
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with  J„  and  real  for  real  parameters.  Thus  a  'n  may  be  determined  by  inspection 
of  an. 

The  coefficients  for  the  radially  symmetric  case  give  g  in  the  form 

co 

(93)  g(r-r0)  =  I  a„cnTn(r-i0),  c„  =  2n  + 1,  2-«5„0,  Tn  =  P„d -fo),  cosn(6 -60). 

n  =0 

The  corresponding  representations  for  <j>  and  u  are 

(94)  </>^l^(kr)incnTn(t-to), 

(95)  u=Z%n(kr)inancnTn, 


and  similarly  the  internal  field  is  the  series 

(96)  <l>=l£n(Kr)inbncnTn. 

The  cases  covered  by  (85)  follow  from  direct  application  of  (4)-(7).  We  may  obtain 
g  of  (93)  from  u  of  (95)  by  using  Wn~i~ntK,  and  we  may  also  obtain  u  of  (95)  by 
substituting  g  of  (93)  into  (25). 

Similarly,  in  order  to  obtain 


(97)  g'(r-fo)  =  X  a'ncnTn(i-i0) 

directly  instead  of  via  the  integral  equation  g'[g]  of  (57),  we  consider  the  analogous 
nonscattering  problems  for  (4)-(7)  with  the  outgoing  wave  form  u  replaced  by  the 
standing  wave  form 

(98)  u'  =  I  LM'ni"a'„cnT„,  Wn  =  i  Im  %n  =  fat?  ~^n\ 

The  asymptotic  form  of  u\  from  ~i~n96  and  dt2)  ~in3€*  is 

(99)  « '  ~  £  a  'ncnTn  (f  •  f0) -k%*Za'ncn(-\)nTn(r- f0)  =  \Xg'(r  ■  f0)  -  \%*g'(-i  •  f0), 

where  we  use  (-l)nrn(jt)  =  Tn(-x)  to  exhibit  the  required  standing  wave  form  (42). 

Corresponding  to  (85)  and  (87)  in  terms  of  (86),  the  internal  coefficients  bn  of 
and  the  analogues  b'„  of  </r',  are  given  by 


(100) 


bn  =  ~iW/$n (Ka )[RnXn  -%'„}  =  b'J(\-a'n), 

b'n  =  -W/Jn(Ka)[RJTn  -rn]  =  -Wl$n{Ka)Dn. 

Using  these  and  an  and  a  we  consider  the  energy  theorems  (34)  and  (60).  From 
oa  =  iA}  =  -Re  \  (il/*dd//ik)  dSP  we  construct 

=  -ST  1 1 bn  1 2c„  Im  [Bvfn  (Ka  )f*  (Ka )] 

=  -yi\bn\lcn\^n(Kaf  Im  Rn, 

where  =  47ra2,  lira,  2,  and  Im  R„  SO.  Similarly  cr'A  corresponds  to  replacing  b„  by 
b'n.  Thus  from  (100), 

(102)  o-k  =  -tfW2 1  cn  Im  RJ\Dn |2  =  -o-0 1  Re  a  '„cn  =  -o-0  Re  g'(k  •  k) 


U01) 


with  StW  =(To,  and  Real,  as  in  (88).  Similarly, 

<ta  =  -yw2 1  cn  Im  RJ\Dn  |2|  1  -  a  ^  |2  =  -mo  I  (Re  a„  + 1  a„  |2)c„ 

(103)  =  -mo  [Re  g (k  •  k) + Jif  |g (f  ■  k)|2] 

=  -<ro  Reg(k'k)-crs. 

in  terms  of  Re  a„  of  (91),  and  <rs  =  |cr0{K*,  m}  =  Re  f  (u*dru/ik)  dif  =  <ToM\g(i -k)|2. 
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As  another  application  of  g[g']  of  (66),  we  generalize  existing  long-wavelength 
approximations  [11],  [12]  for  lossless  scatterers  to  include  loss.  For  lossless  ellipsoids, 
Imgn  to  0(fc5)  was  derived  [12]  by  perturbation  procedures,  and  Reg)2  to  0(k*) 
was  obtained  from  the  scattering  theorem  (32).  Thus  [12] 

g(r2,ii)^ig' ~g*,  g'  -  k3Ai  +  k5A5  +  0(k7),  -g*  =  k6A6  + k*As  + 0(k'\ 

(104)  -A6(r2,  fi)  =  JtA3(i,  r2)A3(r,  f,), 

-Ag(f2,  f i )  r2)A5(r,  r,)  +  As(r,  r2)A3(r,  r,)], 

where  the  A„  were  real  functions  of  real  parameters  B  and  C  =  B172  so  that  g1 
and  gR  corresponded  to  the  imaginary  and  real  parts  of  g.  For  complex  B  and  C,  we 
write 

(105)  g'(?2,ri)  =  ig'(r2,r,) 
and  iterate  (66)  to  obtain 

(106)  g(f2,  fi)  =  g'(?2,  fi)  +^g'(r2,  r)g'(f,  f,)  +  0(g'3), 

where  the  term  0(g'3)  is  0(k9).  Thus,  by  inspection  of  the  existing  results  [12]  for 
form  (104),  we  have  for  complex  B  and  C 

(107)  Im  g(r2,  ri)  =  Re  g',  Re  g(r2,  r,)  =  -Im  g‘  ~gR, 

where  Im  gl  and  gR  correspond  to  absorption  and  scattering  losses  respectively. 
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